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■ Abstract 

This article describes the code SusHi (for "Supersymmetric Higgs" ) [1] which calculates the 
Q \ cross sections pp/pp — > <fi + X in gluon fusion and bottom-quark annihilation in the SM 
and the MSSM, where <fi is any of the neutral Higgs bosons within these models. Apart 
from inclusive cross sections up to NNLO QCD, differential cross sections with respect to the 
Higgs' transverse momentum p T and (pseudo-) rapidity y{rf) can be calculated through NLO 
^ ! QCD. In case of gluon fusion, SusHi contains NLO QCD contributions from the third family 
of quarks and squarks, NNLO corrections due to top-quarks, approximate NNLO corrections 
due to top-squarks, and electro-weak effects. It supports various renormalization schemes 
for the sbottom sector and the bottom Yukawa coupling, as well as resummation effects of 
higher order tan /3-enhanced sbottom contributions. SusHi provides a link to FeynHiggs for 
the calculation of the Higgs masses. 
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1. Introduction 

The recent observation of a new boson at the Large Hadron Collider (LHC) 0, [if has 
opened a new chapter for Higgs phenomenology [if. For the clear identification of this 
particle, precise predictions for Higgs production and decay will be absolutely essential. 
The current status of these efforts is collected in the reports of the LHC Higgs cross section 
working group j^, @] . 

The main production mechanism for a Standard Model (SM) Higgs boson at a hadron 
collider is gluon fusion, where the gluon-Higgs coupling is mediated mostly by virtual top- 
and bottom-quarks (~ 7%). The total inclusive cross section is known through next-to-next- 
to- leading order (NNLO) in quantum chromodynamics (QCD) [7-18]. Even higher order QCD 
effects have been calculated through resummation; electro-weak contributions reach 

up to 8% with respect to the leading order (LO) cross section [25-27 . 

The gluon-fusion mechanism for the neutral Higgs bosons in the Minimal Supersymmetric 
Standard Model (MSSM) is mediated by quarks and their superpartners, the squarks. For 
the CP-even MSSM Higgs bosons h, H, the QCD effects due to quarks can be simply taken 
over from the SM by a rescaling of the cross section with the corresponding modified Yukawa 
couplings. For the CP-odd Higgs boson A, the NNLO QCD corrections to the quark-induced 



total inclusive cross section have been calculated in Refs. [18|, [28|, [29[ . 

Squark contributions to the gluon-Higgs coupling are typically suppressed by powers of 
m q /mq, and thus of importance mostly for small to moderate squark masses m,. For the 
CP-odd Higgs boson A, they are absent at LO0 



Comprehensive reviews of the Higgs theory within the SM and the MSSM can be found in Refs. 
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The individual components for a calculation of NLO QCD corrections in the MSSM are 
known: The NLO cross section due to (s)top (s) quarks, gluons, and gluinos was calculated 
for a CP-even Higgs boson with mass in Refs. [32-34 ] by applying an effective-theory 
approach in the limit <C m t , mi, m g , similar to what is used in the SM at higher orders. 
In this approach, even NNLO effects have been first approximated 35], and recently fully 



calculated [36|, |37|. The analogous NLO result for a CP-odd Higgs was first obtained in 
Ref . F 



. Due to the smallness of the bottom-quark mass, these results cannot be transferred 
to the bottom/sbottom sector. However, the limit m^, m;, <C mi, m g is applicable in a large 



region of the parameter space, and the corresponding results were presented in Refs. 39 



40 



Recently, more general results for the cross sections, allowing for larger Higgs masses, 
were obtained for the CP-even and -odd Higgs bosons and both for the top/stop and the 
bottom/sbottom sector, in Refs. [4l|, |42|. A fully numerical calculation of the gluon- Higgs 
form factor for general quark/squark/gluino/Higgs-masses has been reported in Ref. 43 . 



For the pure squark contributions, the full Higgs-mass dependence of the NLO contribution 
to the cross section was presented (numerically and/or analytically) in Refs. 44|-|46 . 

For large values of the MSSM parameter tan /3, the coupling of the light CP-even Higgs 
boson to bottom-quarks is significantly enhanced relative to the SM Yukawa coupling, so 
the bottom sector may be much more important for gluon fusion in the MSSM. In addition, 
an enhanced bottom-Higgs coupling increases the cross section of another Higgs production 
mechanism in the MSSM, namely associated production with bottom-quarks, pp/pp — > bb<f). 
If the final state quarks are not tagged, a suitable theoretical approach to the cross section 
is the process bb — > <fi, called bottom-quark annihilation in what follows!! It resums terms of 
the form lnmj/m^ by means of 6-parton distribution functions (PDFs) and was calculated 



up to NNLO QCD in the SM [48|, |49j. The result can be directly translated into the MSSM 



by rescaling it with the proper bottom- Yukawa couplin g; e ven the dominant sbottom effects 



50 



51 . SusHi evaluates both the cross 



can be taken into account by an effective coupling 
section for gluon fusion and bottom-quark annihilation. 

For gluon fusion, SusHi includes results for all NLO QCD contributions due to the third 
generation of quarks and squarks. The real corrections at NLO are well-known; SusHi 
implements them using the routines of Ref. [40j. For the virtual corrections to the pure 



quark diagrams, it uses the analytic expression of Ref. 52) which was obtained from the 
integral representation in Ref. 10|. Concer ning the genuine virtual supersymmetric (SUSY) 



corrections, it employs the results of Refs. |32|, |33j, |41, 
and the (s)bottom-mediated gluon-Higgs coupling, respectively. 



42| and Refs. |39j, |4jJ for the (s)top 



NNLO QCD__effects are 
3, and 
531. Electro- 



11 13, 28 



taken into account for the top-quark induced gluon-Higgs coupling 
approximately for the top/stop/gluino- induced one [35j by using gghOnnlo 
weak corrections 25-27[ are included as tabulated correction factors. The cross section is 
provided in various renormalization schemes (in particular in the sbottom sector), allowing 
for an on-shell, DR, or a dependent renormalization of the soft-breaking parameter A^, for 
example. In addition, the bottom- Yukawa coupling can be chosen on-shell or in the MS- 
scheme. Higher-order sbottom effects can be included through the parameter 54]-|59 . 



2 See Ref. [47[ for a more detailed discussion. 

3 



For the calculation of the bottom-quark annihilation cross section, SusHi makes use 



of bbh@nnlo [60] and re-weights its results by the MSSM couplings. It uses the LHAPDF 
library |6l| which allows to conveniently switch between different PDF sets, and it can be 
linked to FeynHiggs 62h65| for the two-loop calculation of the Higgs boson masses in the 
MSSM. 

Apart from inclusive cross sections for gluon fusion and bottom-quark annihilation, SusHi 
allows for (upper and lower) cuts on the transverse momentum and/or the (pseudo-) rapidity 
of the outgoing scalar cf). In case of gluon fusion, differential distributions with respect to 
these kinematic variables can be obtained (for not too small). 

Note that a number of codes for the calculation of Higgs cross sections in the SM and 
the MSSM exist, see Refs. 66M72l| . for example. They overlap with SusHi to a greater or 
lesser extent; the distinctive feature of SusHi is to provide full NLO QCD (and partial 
NNLO and electro-weak) corrections for the dominant production mechanisms of the three 
neutral Higgs bosons of the MSSM, both inclusive and differential, in various renormalization 
schemes. Further details will be given below. 

The remainder of this paper is organized as follows: In Section [2J we present the physical 
background of SusHi, recalling the framework of the Higgs and quark/squark sectors in the 
MSSM with special emphasis on the renormalization of the squark sectors and the resumma- 
tion of tan /^-enhanced sbottom corrections in the bottom Yukawa coupling. Subsequently, 
we discuss the various contributions for the calculation of the gluon-fusion and the bottom- 
quark annihilation cross section as they enter in SusHi. We briefly describe the kinematic 
variables for which cuts can be applied and distributions be obtained. In Section El we 
describe the program SusHi, in particular its workflow, installation, and usage, as well as 
the input and output files. Our conclusions are given in Section [71 Appendix A contains 
the couplings of the squarks to the Higgs bosons 0. 



2. Physics background 

This section first introduces our notation for the relevant parts of the SM and the MSSM. 
It describes the renormalization of the squark sector and the possible choices for the bottom 
Yukawa coupling provided in SusHi. 

2.1. Standard Model 

The SM contains one scalar weak isospin doublet which corresponds to a single physical 
particle H (electric charge and spin zero). The pure Higgs sector of the Lagrangian is 
determined by two parameters: the vacuum expectation value v ~ 246 GeV, and the mass 
of the Higgs boson rajj. The Yukawa couplings of the fermions to the Higgs boson are given 
by Yf = \/2mf/v, where my is the fermion mass. SusHi requires all input masses in the 
on-shell scheme, except for the charm- and bottom-quark mass which has to be given as 
the MS mass m q (m q ) = (to^ ), q e {c, b}. The calculation of various internal bottom 
masses is addressed in Section 12.31 
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2.2. Supersymmetry 

The MSSM contains two Higgs doublets, named and H u , which develop the vacuum 
expectation values = v cos/3 and v u = v sin/3, where the parameter /3 is undetermined. 
They form two CP-even Higgs fields h,H, one CP-odd (or "pseudo-scalar") Higgs field A, 
and two charged Higgs fields H ± . At lowest order, the mass spectrum of the Higgs sector is 
determined by SM parameters, tan/3 = v u /vd, and the CP-odd Higgs mass vtia- Radiative 



corrections to the Higgs mass spectrum are generally quite large |73l-l75| ; currently, they 



are known through three-loop order |76H78| . SusHi provides a way to conveniently take 
into account two-loop corrections by linking to the program FeynHiggs [62 65]. This is an 
optional feature, however; any Higgs mass can be given as an input to SusHi. 

The mixing of the isospin to the Higgs mass eigenstates in the CP-even sector is governed 
by the angle a. Together with /3, it determines the relative strengths of the Higgs boson 
couplings (ff (</> G {h, H, A}) to the SM fermions with respect to the SM Higgs boson cou- 
Appendix A ), which thus enter the fermion Yukawa couplings Yjt = \/2mfg^/v: 



plings 



sec 



if u 



9 h d 



cos a 
sin /3 
sin a 



9u 



9d 



sin a 

sin /3 ' 
cos a 



at 



9d 



tan/3 ' 
tan P . 



(1) 



cos/3 "" cos/3 

These normalized couplings are independent of the fermion generation. Our calculation 
includes the third generation of squarks which enters the Lagrangian in the form 



with the mass matrix 



Mi 



M 2 +m 2 q + m| cos(2/3)(T 9 3 - Q q s 2 w ) 

m q (A q - flK q ) 



m q (A q - fjLK q ) 
M$r + m 2 q + m 2 z cos(2/3)g 9 s^ 



(2) 



(3) 



for an arbitrary species of squarks q. This formula contains the SUSY soft-breaking parame- 
ters Ml, M~ R , and A q . The parameter \i determines the mass of the fermionic Higgs partners, 
the Higgsinos, whereas the Z-boson mass mz and the weak mixing angle 9w (sw = sin^v^) 
are the usual SM parameters; m q , Q q , and T q are the mass, the electric, and the weak charge 
of the corresponding quark q, respectively. It is K h = tan/3 =: tp and n t = 1/tp. 

The physical particle states are obtained by the diagonalization of the mass matrix in 



Eq. ([3]) which we do in accordance with Ref. 79| using 



with 



U~„ 



cos 9 q sin 9 q 
— sin 9 r -, cos 9; 



(4) 



By choosing < 9 q < ir, the masses of the squarks q~i and q 2 are ordered m ql < m q2 and 
given by the square roots of the eigenvalues of .M- in Eq. (J3J): 



m 



ql2 



(M 2 L + M~ R 

1 
2 



+ m 2 q + -T^m 2 z cos(2f3) 



(5) 



T -^(M 2 + M 2 R + m% cos(2/3)(T3 - 2Q q sin 2 9 W )) 2 + Am 2 q {A q 
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jlK q 



The entries of A4| can also be expressed in terms of the mass eigenvalues and the mixing 
angle 

2 = (ml cos 2 6 9 + m% sin 2 9~ q (mf, - m 2 2 ) sin 9 9 cos 

* V( m |i — m \2) sin^cos^ mpsin 2 ^ + ni~ 2 cos 2 9 q J ' 

which implies that the mixing angle can be obtained from 

sin(20 3 ~) = 2mg( 2 ^~f g) . (7) 

Note that Eq. (J7J) does not uniquely define 9 q ; a shift #q — )■ | — 9 q , which corresponds to 
sin^g cos9 q , might be in order to allow for m q \ < m q2 . For completeness, the couplings 



of the squarks to the MSSM Higgs bosons can be found in |Appendix A 



2.2.1. Renormalization of the (s)top sector 

Using Eq. (J2J), we eliminate A t from the (s)top sector contribution of the amplitude before 
renormalization and express it in terms of the on-shell parameters for the top mass mt, the 



stop masses and mf 2 and the mixing angle 0j, defined according to Section 3.1 in Ref. [79 
In practice, the user specifies the soft-breaking parameters Ml = Mi(t), M^ R , and A t , as 
well as the on-shell top-quark mass m t ° . Setting m t = m t os , SusHi inserts them into the 
mass matrix of Eq. (EJ) whose eigenvalues m 2 v m~ 2 , as well as the corresponding stop mixing 
angle 9 t are interpreted as on-shell parameters^ Note that in case FeynHiggs is used, the 
on-shell stop masses and the stop mixing angle are simply taken over from its output. 

2.2.2. Renormalization of the (s)bottom sector 

The renormalization of the (s)bottom sector is more subtle. At tree-level, the soft- 
breaking parameter M 2 is identical for the sbottom and stop sector due to SU (2) L symmetry. 
At higher orders, however, in the on-shell scheme we distinguish M 2 (t) and M 2 (b) similar 



to Ref. |80H83| by 

M 2 (b) = M 2 L (t) + 5M 2 (t) - 6Ml(b) = M 2 (t) + AM 2 , (8) 
with the individual counterterms given by 

5M 2 L {q) = cos 2 9 q 5m%f S + sin 2 9 q 5m 2 - 2 os - {m\ x - m 2 2 ) sin(2# 9 ~)<^? s - 2m q 5m° s , (9) 

where the on-shell counterterms 5m^ os , <5m~ 2 os , <5m° s , and 8 99 s are analogously defined as 
in the top sector, see Ref. [79 . 

Note that the finite shift AM| depends on m~ hl ,m~ b2 , and 9b- In order to determine its 
numerical value, we first calculate "tree-level" values for the sbottom masses and mixing 
angle by inserting the parameters 

Ml = (M2(6)r e = M 2 (t) , M 2 R , A b , m° s (10) 



3 This can be seen as an indirect definition of the renormalization scheme for Ml, Mf R , and A t . 
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into the mass matrix ([3]) (setting m& = m® s ). All parameters of Eq. (11 01) are input to 
SusHi, except for the on-shell bottom-quark mass m° s which is determined from the input 
parameter mb(mb) as described in Section 12.31 With these tree- level sbottom masses we 
define the scale 

1 



/id 



:(m~ g + m~ bl + m~ b2 ) 
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As was pointed out in Refs. 39, 7§, 84, 85 1, replacing A b in the amplitude through Eq. (J7J) 
before renormalization - analogous to the stop sector - leads to potentially large corrections 
SAf, cx (a s /i 2 ta.n 2 /3) /rrig. It was therefore suggested to use Eq. ([7]) in order to eliminate nib- 
The counterterm for rrib in this scheme, denoted 5m b ep , is obtained from Eq. ((7|): 



Sm b ep = 2m b cot(29 b )S9 b 



2m 2 b ■ SA b 



sin(26b)(m 



m 



+ m b - 



r 2, OS r 2,OS 

om- — om 7 



bi 



62 



'12) 



i/i > 



m 



61 



m 



62 



Here it is already implied that we always renormalize the sbottom masses on-shell, while 
the renormalization of Ab and 9b is still unspecified. 

In order to calculate the on-shell sbottom masses, we choose the "on-shell" renormaliza- 



tion of Ab, defined through a kinematical condition on the A&i^-vertex |39l.l79l. l84l. l85j|. The 
corresponding counterterm is 



5A% S = (A b + li cot /3) 



X dep 

ft 2 2 \ i el 2 2 \ 0m b 

f(m- bV m- b2 ) + f(m- b2 , - — — 

I lib 



(13) 



with 



f{ml,ml 



q 8 (^d) 2 

3 



7T 



4 + 2 log 



/4 



mt 



m g 



BQ n (m(,m b ,mg, fid) 



A h + jx cot /9 



:i 



m 2 



mt 



m\ 



(14) 



where the function B^ n can be taken from Eq. (B.8) in Ref. [33] replacing //r — )■ /i^. Solving 
Eq. (El (with cL4 b = <5A° S ) and Eq. flT3J for 5A° S , we obtain 



8Af s = (Asin/3 + /icos/3) {-(m^ - mf 2 ) sin(20 6 ) [/(mgj, mf 2 ) + /(mf 2 , mfj] 



+ (5m| os - 5m- 2 2 OS ) sin(20 fe ) + 2^ 6 (mf 1 



62 



cos(2# fe 



(15) 



[2 (sin /3 (A 



bm b 



m 2 — m- b „) sin 6*6 cos #&) + cos 



-i 



which in turn allows to calculate 5m b ep . Note that one is still free to choose the renormal- 
ization condition for 9 b - 

The on-shell sbottom masses are finally obtained as follows: We calculate the counterterm 
Sm b ep from Eq. f JT2|) using SA b = 5A% S and 59 b = 50? s , the "tree-level" sbottom masses and 
mixing angle and the on-shell bottom mass. This yields a numerical value for AM| and 
similarly for m b ep = mP R (/id) — 5m b cp + 5m^ R , where the DR mass is calculated in the SM. 



Diagonalizing the mass matrix .M- in Eq. fl3]) with m& = m b° P i we thus obtain on-shell 
sbottom masses m bl and m~ b2 and the on-shell sbottom mixing angle 9~ b . In case FeynHiggs 
is called for the calculation of the MSSM Higgs masses, we take m b ep as well as AM| from 
this program which, when inserted into the sbottom mass matrix Eq. (J3]), results in on-shell 
sbottom masses and mixing angle consistent with the values given by FeynHiggs itself. 

After the on-shell sbottom masses have been determined, SusHi allows for a change 
between various schemes regarding the renormalization of mj, and 9~ b in the sbottom 
contribution to the gluon-fusion amplitude. Note, however, that either mj or Ab is required 
to be a dependent quantity; a dependent renormalization of the sbottom mixing angle is not 
offered as option. The possible choices are summarized in Tab.[TJ Numerical differences in 
the cross sections between the various renormalization schemes and their implications will 
be investigated in a separate publication. 



Scheme choices 


m b 


dep. 


OS 


DR 




dep. 


OS 


DR 


H 




OS 


DR 



Table 1: Available renormalization schemes for the sbottom sector. The default option is marked in red. 

To summarize, apart from fi, tan /3 and the SM parameters, the input parameters of 
SusHi that determine the squark sectors are: 

• M\, M~ R , A t , which directly determine the on-shell stop masses and mixing angle 
through diagonalization of Eq. ([3]) 

• M~ R , Ab, rrib(mb), where Ab is understood as renormalized according to Eq. ([TBI . 

Switching between renormalization schemes changes both the counterterms to the amplitude 
as well as the numerical values of the parameters nib, m ii^ m &2' Ab, and 6 b , of course. We 
remark that changing the renormalization scheme in the sbottom sector affects nib only 
in the Higgs-sbottom couplings; the renormalization of the bottom mass occurring in the 
Higgs-bottom Yukawa coupling is independent of that (see Section [275]) . For the bottom 
mass occurring in internal propagators (rather than in couplings), SusHi always uses the 
on-shell value. 

The switching between the different renormalization schemes is done at the renormal- 
ization scale /ir to guarantee the same coupling strength a s (/iR,) for the NLO counterterms 



and the cross section itself. In case SusHi makes use of the formulas in Refs. [3 91 . |41[, the 
counterterms at NLO are expanded to the correct order in the bottom mass to match the 
expanded NLO amplitudes. 
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2.3. Bottom mass calculation /Resummation of tan f3- enhanced contributions 

The bottom-quark mass in the input files of SusHi is inserted in the MS scheme m b (m b ). 
Together with the input value of a. q (mz) : we calculate ot s {m b (m b )) by 4-loop running with 



5 active flavors. Using Eq. (13) of Ref. [86J] at 3-loop level (see also Ref. [87(] ) , m b (m b ) is 
transformed into its on-shell value m^ s . 

In the SM, the bottom Yukawa coupling Y b = y/2m^ jv can be chosen m% = m® s or 
alternatively m( = m^ s (/i 6 ), where \i b G {m b ,fi R }; /x R denotes the renormalization scale. 
As indicated above, the bottom mass entering the internal propagators is always set to the 
on-shell mass m° s in SusHi. In the MSSM, SusHi offers various options regarding the choice 
of the Higgs-bottom Yukawa coupling = y/2n\'^ / v : 

on-shell coupling: m b = m^ s (16a) 



basic resummation: mY'^ = b —— (16b) 

b 1 + A b ( J 

full resummation h: mY ,h = b —— ( 1 — A b ] (16c) 

b 1 + A b V tan/3 J v ; 

„ Y,H ™? S f, . a tan oY 

m h ' = - — 1 + A b (16d) 

b 1 + A b V tan/?/ v ; 

y,a_ m b S ^ - 1 



OS 



A ' b = T+Al V 1 ~ A *tan^8 ) < 16e) 

Running coupling: mj'^ = ^ — (16f ) 

1 + A b 




Herein, A& resums higher order sbottom contri- 
butions as shown in Fig.[T] [54j-l59| . for example, 
calculated at the scale /id, defined in Eq. fTTT]) . 
The exact formula is given by: 



A, 



Figure 1: Feynman diagram inducing an effective 

coupling of the bottom-quarks to H u . The coupling /(a, 6, c) 

is proportional to [itp and can be resummed in A5. 



—a s (/i A )mg/itpI(m 



2 2 2\ 

~ bV rn i2 ,m~ g ) 



a61n(f) + 6cln(|) + caln(|) 
(a — b)(b — c)(a — c) 



(17) 



If the input values for SusHi are determined by FeynHiggs, SusHi allows to use the value 
of A b as given in the output of this program which contains also the electro-weak contribu- 
tions from neutralinos and charginos in accordance with Ref. 88]; two- loop corrections to 
Ab |89M91| are not yet included. For the running coupling in Eq. (116fj) . the scale \i b can be 



set to m b or /iR. The numerical differences between the various schemes will be discussed in 
a forthcoming publication. 
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3. Cross section for gluon fusion 



As indicated in the Introduction, the most important production channel in the SM and 
for moderate values of tp in the MSSM is gluon fusion. After quoting the well-known results 
for the LO cross section, our implementation of the NLO contributions is explained. A 
discussion of NNLO and electro-weak contributions follows. 

Using the notation of Ref. [10j, the hadronic cross section for <\> e {h, H, A} at NLO QCD 
can be written as follows: 

a{pp ^ (j) + x)=at[l + C+^\ t^— + Aa* g + Aa* q + A4 , (19) 

where = m^/ s, with the hadronic center-of-mass energy s. The factor ctq is determined 
by the LO cross section, arises from NLO terms in the partonic cross section that are 
singular as s — > (s is the partonic center-of-mass energy), and 

-g(x)g(r/x) (20) 



dr J T x 

is the gluon-gluon luminosity. The quantities Acr^ 9 , Acr^, and Acx^ comprise the terms that 
are regular as s — > m| in the partonic cross section, arising from gg, gq and qq scattering, 

respectively. Loosely speaking, is due to the virtual, while the Acr^ are due to the real 
radiation contributions. The latter are implemented in SusHi by expressing them in terms of 



Passarino-Veltman functions [92[|, see Ref. 40]. The strong coupling a s for the cross section 



calculations is renormalized in standard QCD with five active quark flavors. 
3.1. LO cross section 

Taking into account the third generation of quarks (and squarks in the MSSM), the 
normalization factor in Eq. (I19p can be written in the form 



\A^ 0) \\ (21) 



" 288^ 

with Fermi's constant Gf- The amplitude A for is given by 

A+>® = « (0) + < {0) ) , (22) 

q£{t,b} 

with the individual contributions for 6 {h, H} 

< (0) = 0^(1 + (1 - rt)f{rt)l < (0) = ~ f ^4,(1 - r*f(4)) , (23) 

8=1 

using the notation 

g ~ my v ~ ml [ } 
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and the function 



arcsm 



2 1 



f(r) 



in 



T > 1 
T < 1 



For the CP-odd Higgs = A, the squarks do not contribute at LO, i.e. 
quark contribution can be written in the form 



(25) 



0, and the 



,A(o) 



3r A 



(26) 



The couplings of the Higgs <\> to the quarks can be taken from Eq. (jTJ), the couplings g'-^ 
to the squarks from Appendix A Needless to say, in the SM, the squark couplings have to 
be set to zero and the quark couplings to g% = 1. 

3.2. NLO virtual contributions 

As indicated above, the coefficient contains the virtual corrections to the gg process 
and is regularized by the infrared singular part; moreover, it includes the counterterms to 
LO quantities. We write it as 



C 4, = 2Re 



^,(i) 



,(o) 



+ vr 2 + /3 log ^ , (27) 

where /3q = 11/2 — n//3 with rif = 5; /xp and /xr, denote the factorization and the renormal- 

3ctively. The 1 
limit of large stop and sbottom masses are given by 



ization scale, respectively. The NLO amplitude A?'^ and the LO amplitude Ato in the 



Mi) 



A<PM 



q&{t,b} 




(28) 



Available results for the NLO contributions have been discussed in the Introduction. In 
SusHi, we use the analytic formulas of Ref. 52] for the quark-induced terms af . The purely 



squark-induced terms (see Fig.|2](a), for example) need to be considered in combination with 
the mixed quark/ squark/gluino diagrams (two examples are shown in Fig. [2]) in order to 



preserve supersymmetry, resulting in the coefficient Oq 
these amplitudes in two limits: 



SusHi implements expansions for 



<C m q , rriqi, mp, m g [321434] which is valid through < min(2m g , 2mg, rrig + m q + 



and thus applies to the top-stop sector as long as is not too heavy. SusHi 
incorporates the publicly available program evalcsusy.f [93J in order to use this 
result for the light Higgs h. 



m^m q <C m q x^m q 2,m,g [39r|42| which holds through < mm(2mq,m q + trig), and 



thus applies to the bottom-sbottom sector, as well as to the top-stop sector when 
heavy; SusHi uses the formulas of Refs. 4l|, |42| in these cases. 

11 



is 



Note that both expansions hold only as long as the Higgs mass is not much heavier than 
the typical SUSY mass. For larger Higgs masses, only the fully numerical result of Ref. 43[ 
is known so far. 




(a) (b) i (c) 

Figure 2: Example Feynman diagrams contributing to : gluon-squark contribution (left), gluino-quark- 
squark contribution (middle) and gluino-quark-squark contribution (right). The latter is partially resummed 
by the usage of Af, from Eq. (p~7)) . 



3.3. NNLO corrections 

While the implementation of NLO corrections allows for the evaluation of inclusive and 
exclusive quantities, NNLO corrections are available in SusHi only for the total inclusive 
cross section. In addition, only the top- (s) quark induced gluon-Higgs coupling is taken into 
account at NNLO; the top-squark induced one only approximately: 

MSSM _ MSSM , I t _ t \ (r,Q\ 

u gg4>,NNLO ~ u gg<p,NLO ' \ u gg0,NNLO u gg<f>,NLO) J ) 

where, on the right-hand side, the NNLO term is evaluated with the NNLO PDFs, while the 
NLO terms are evaluated as usual with NLO PDFs. 

The (N)NLO top-(s)quark contributions to the cross sections <r*^ (n)nlo are calculated 
with the help of the programs ggh@nnlo [53] and evalcsusy [93] which work in the effective 
theory approach of heavy top (s) quarks. The NNLO top-squark (and mixed top/stop/gluino) 
effects, available for (f> = h only, have been evaluated 36, 3?| by applying the limit <^ 



m q ,mqi,rriq2,mg. The result has been implemented in a computer code that involves 
Mathematica and a number of other programs. SusHi includes an approximation of these 
NNLO effects according to Ref. [35j0. 

3.4- Electro-weak corrections 



The full NLO electro- weak (EW) corrections are known only in the SM 25|. It has been 



suggested to assume complete factorization of QCD and EW effects |94|, thus writing 



a 



SM,t , x ^SM,t 

ggH,NNLO,EW " 



1 + ^Ew) Cr < ; < ; i/ i NNLO ' (30) 



4 The NNLO term of the normalized Wilson coefficient (k2 in the notation of Ref. [35]) is replaced by its 
SM value in SusHi. 
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For the CP-even Higgs bosons in the MSSM, a formula based on Eq. (|29|) and Eq. ( 130]) 
has been used for the combination of QCD and electro- weak corrections in Ref. 401 ] : 



a 



MSSM 



gg</>,NNLO,EW — a gg4>,NLO 



MSSM 



+ (1 + 5 EW )a 



EW ) U 33</>> nnl O °"ss0,NLO J 



or, at NLO precision, 



MSSM 

'gg<j),NLO,EW ~ u gg<j),NLO 



, r • . MSSM _ r t 



(31) 



(32) 



Alternatively, it has been suggested in Ref. [67| to use the SM electro-weak corrections 
due to light quarks only 26, 27]. Following Ref. [671 ]. we define the correction factor 



°EW 



a EM 2Re (A^A*' EW ) 



71 



(33) 



where A^'^ denotes the complete LO amplitude including quark and squark diagrams, see 
Eq. ( 12"2"|) . and the electro- weak amplitude is given by [27 



A<j>,EW _ ^ Xw 

~ 8 s 2 



2 /5 7 22 



with 



Xy 



m 



1 / r y 

2 , my — z — 



Ve{W,Z}, 



(34) 



(35) 



the electro-magnetic coupling oem, and sw = sin6*iy = (1 — c 2 ^) 1 / 2 the sine of the weak 
mixing angle. Supersymmetry enters through the relative couplings g v given by 



g v = sin(/3 - a) 



gy = cos(/3 - a) 



(36) 



The function A\ [x] can be found in Refs. [26|, |27J . Since its numerical evaluation is rather 
involved, SusHi implements (and S EW ) in terms of an interpolation grid in (and 
m t ), using fixed values for the gauge boson masses and widths, as well as for the weak 
mixing angle (m w = 80.385 GeV, T w = 2.085 GeV, sin 2 6 W = 0.22295, m z = 91.1876 GeV, 
Y z = 2.4952 GeV j95|); the input values to SusHi for these parameters are ignored in the 
evaluation of the electro-weak corrections. The electro-weak correction factor due to light 
quarks 5{f w multiplies the NLO MSSM cross section, while the NNLO QCD effects are simply 
added as in Eq. (129"]) : 



a 



MSSM 
99<t> 



MSSM /i , Hf \ , t 
a ggcj>,NLO\ l °Ew) ~t~ ° g 



sNNLO 



sNLO 



(37) 



SusHi leaves it up to the user to decide whether to use Eq. (|3"T|) or Eq. f l3T|) in order to include 
the electro- weak corrections. For a SM-like Higgs and < 2m t , both approaches lead to 
comparable NLO results. The EW corrections for a CP-odd Higgs are not known and thus 
not included. 
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4. Cross section for bottom-quark annihilation 



In supersymmetric theories, where the Higgs coupling to bottom-quarks can be enhanced 
by tan/3, associated production (bb)4> + X can be similarly or even more important than 
gluon fusion. Two theoretical approaches have been pursued for the theoretical description 
of this process: In the four-flavour scheme (4FS), the relevant production processes at lowest 
order QCD are gg — > (bb)<fi (see Fig.[3](a)) and quark-antiquark annihilation qq — > {bb)<p (96^ 



98] . However, when integrating over all final-state bottom-quark momenta, potentially large 



logarithms lnm^/m^ occur. They can be resummed b y th e introduction of bottom-quark 



PDFs, which defines the five-flavour scheme (5FS) [99l. 1 1 Of ] . The LO process in this latter 
scheme is bottom-quark annihilation 66—7-0 for which the lowest-order Feynman diagram 
is shown in Fig.|3](b). 




(a) (b) 

Figure 3: Feynman diagrams showing the associated production process gg 
bottom-quark annihilation in the 5FS (right). 



(bb)(f> in the 4FS (left) and 



SusHi implements results for associated 



section, it links the program bbh@nnlo [49 



-production in the 5FS. For the inclusive cross 
in order to obtain the NNLO QCD prediction 
for the bottom Yukawa c oupl ing; . This is then re-weighted by the 



corresponding resummed SUSY coupling g£ [58|, 1101] as follows 



a 



MSSM 



9? 



SM 
bbH 

9? 



(9. 



b) 



1 + A fc 



1 + A, 



with 



tana 
'tan/3 



~9 h b 



9 h b 



1 + A 6 



9b 



9t 



cot a 



tan (3 



1 + A 6 



1- A, 



tan 2 (3 



(38) 
(39) 



where the g$ are given in Eq. ([TJ and A b is determined as described in Section | 

For differential cross sections due to bottom-quark annihilation, SusHi includes the NLO 
virtual corrections for bb — > <p and combines them with the LO real-radiation processes 
bb — > g<p and bg — > b(j) using dipole subtraction!! Similar to the fully inclusive case, we 
multiply with the resummed SUSY couplings to obtain MSSM cross sections. 



5 We are grateful to M. Wieseman n for pro viding us with the corresponding Fortran routines which 
entered the studies presented in Refs. 102i . Il03| . 
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5. Differential cross sections 



Apart from the total inclusive cross sections due to gluon fusion and bottom-quark 
annihilation, SusHi also allows for the computation of differential cross sections in these 
processes. In particular, one may apply upper and lower cuts on the Higgs' transverse 
momentum pt, its rapidity y or its pseudo-rapidity 77, where 

, ( 8\ 1 { \p\+Pl \ 1 ( E + p L \ .... 

Here, p = Pt+Pl is the 3-momentum of the Higgs boson, p^ the longitudinal component, E 
the Higgs boson's energy, and 9 the scattering angle (all in the hadronic reference frame). For 
gluon fusion, SusHi also provides the differential quantities da/dpr, da/dy, and d 2 a/ (dprdy) 
(or, alternatively, drj instead of dy). We add that, since the distribution in y and t] is 
symmetric, minimal and maximal values for y are understood as < y min < \y\ < y max (and 
similarly for if). In order to get reliable results, the precision for the numerical integration 
in SusHi should be set to a higher value for differential quantities than for inclusive cross 
sections. 

Note that at LO, i.e., 0(a 2 s ) for gluon fusion and O(a s ) for bottom-quark annihilation, 
the Higgs' transverse momentum is always pr = 0. SusHi provides results for non-inclusive 
quantities through NLO, i.e., 0(a^) for gluon fusion and 0(a s ) for bottom-quark annihila- 
tion. Let us also add that p^-cuts or p^-distributions should not be too low (pT/ m </> ^ 0.1), 
since otherwise potentially large logarithms may spoil the perturbative convergence of the 
fixed-order results i mplemente d in SusHi. For the resummation of such terms in Higgs 



production, see Refs. [1041 . 1 1051 ] . for example 



6. The program SusHi 

This section describes the most important technical details of the program SusHi, in- 
cluding its installation and usage. 

6.1. Workflow 

The workflow of SusHi is de picted in F ig;.[H The input is controlled by a single input file 



whose format is SLHA-inspired [1061 . 1107] . In case of the MSSM, the user specifies whether 
the Higgs mass is calculated by FeynHiggs or provided by the user himself. After the 
initialization of internal parameters which are derived from the input data, SusHi transforms 
them to the specified renormalization scheme and determines the resummation of tan /3- 
enhanced terms in the bottom Yukawa coupling, see Sections 12.2.21 and 12. 31 Afterwards, the 
gluon fusion and bottom-quark annihilation cross sections are calculated up to the desired 
perturbative order. The NNLO top- (s) quark induced and the electro- weak contributions for 
gluon fusion are taken into account only for the inclusive cross section. Not shown in the 
workflow is the link to LHAPDF which occurs at various stages of the internal calculation. 
The output is printed to the screen and written to an output file which follows the same 
format as the input file. Details concerning the in- and output files are given in Section 16.41 
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SusHi input file in SLHA-style 



SM, SUSY: Higgs 
mass as input 



SusHi initialization 



SUSY: Higgs mass 
by FeynHiggs 



Transformation to specified renormalization scheme; calculation 
of couplings (including resummation) 



Total cross section 



7 f 



bbhOnnlo, bbh diff. 



Diff. cross section 



AT 



Calculation of ggh/bbh at (NN)LO 



Calculation of ggh/bbh at (N)LO 



ggh@nnlo, electro-weak contr. 



Re-weighting ggh: gghOnnlo top 
contribution, electro-weak contr. 



Performing cuts in p?, y, r) 
ggh: da/dp T /dy/di] 



Screen output /output file in SLHA-style 



7 



Figure 4: Internal workflow of SusHi. Red boxes indicate interaction with the user, who has to provide an 
input and gets an output file, if no error messages are shown. Green boxes refer to external code (see text), 
which is linked to/included in SusHi. 



6.2. External code 

As already mentioned, SusHi includes existing code like gghOnnlo, bbh@nnlo, or the 
electro-weak grid. The integration of these programs does not require any action from 
the user though; they are part of the distribution and are simply linked to SusHi upon 
compilation. 

However, SusHi can/must be linked to the following external code: 

• FeynHiggs [62M60I] : For the calculation of the supersymmetric Higgs masses G 
{h,H,A}, SusHi can be linked to FeynHiggs. Its input is controlled via the SusHi 
input files. Note that the current version of SusHi does not support complex MSSM 
parameters, so FeynHiggs is called with the default flags for the real MSSM. 



LHAPDF 



61] : SusHi has to be linked to LHAPDF which provides a large variety of 
different PDF sets. This allows one to change the PDF set used by SusHi simply by 
changing the input file. 
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6.3. Installation and Usage 

A tarball with the source files of SusHi can be obtained from Ref . [lj . Unpacking results 
in a main folder with the following subfolders: 

bin : contains the executable program sushi after compilation 

example : various example input files to be used in bin 

include, lib : object files and libraries 

src : SusHi source files, including external code 

The file README in the main directory contains installation instructions and the history of 
the code. Compilation is most easily done by adjusting and running the Makefile in the 
main folder as follows: 

• Specify the location and the name of the library of LHAPDF, for example: 
PDFLIBP = /usr/local/lib 

PDFLIB = -1LHAPDF 

• If FeynHiggs should be used, specify the main directory of the compiled FeynHiggs 
library, for example: 

FHPATH = /home/. . ./FeynHiggs-x.x.x 

Note that SusHi requires FeynHiggs version 2.9 or higher. 

• Run the configure script in the main folder: 



. /configure 



It tries to find the gfortran or ifort compiler and the dependencies on your local 
machine. If you prefer a different compiler, or if the script fails, you can specify the 
relevant variables (F77 and LDFLAGS) in the file compilerissues yourself. 

In the main folder, run make: 



make [option] 



This command takes either of two optional arguments: make predef =NoFeynHiggs 
omits the link to FeynHiggs which cannot be used in this case; make clean deletes 
all object files, libraries, and the executable sushi. 

After compilation, one may perform a test run of SusHi by copying one of the input files 
from the example- to the bin-folder, change to the bin-folder, and run 



/sushi sushi . in sushi. out 



Note that apart from the input filename, the user also has to provide a name for the output 
file. For parameter scans , we recommend the auxiliary routines SLHAroutines, which can 
be downloaded from Ref. [1081 ] . In the following we will discuss the input and output files in 
more detail. 
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6-4- Input and output files 

SusHi allows for calculations in the SM as well as in the MSSM. Although the former is 
in many respects a limiting case of the latter, SusHi distinguishes both cases, and we will 
discuss them separately in what follows. 



6.4-1- Standard Model 

A typical input file for a SM calculation is shown in the following. It divides into blocks, 
each of which contains a number of entries, specified by one or more leading blanks, an 
integer "entry number" , a value, and a comment initiated by the hash symbol #. 



Block 
1 
2 
3 
4 
5 
6 
7 

Block 
1 
2 
3 
4 
5 
6 
8 

Block 
1 

Block 
1 



SUSHI 




8000. dO 

2 

2 

1 

SMINPUTS 
1 



# model : = SM , 1 = MSSM 

# = scalar Hlggs (h) , 1 = pseudoscalar Hlggs 

# collider: = p-p, 1 = p-pbar 

# center -of -mass energy in GeV 

# order ggh : -1 = of f , = L0 , 1 = NLO , 2 = NNLO 

# order bbh : -1 = of f , = L0 , 1 = NLO , 2 = NNLO 

# electroweak cont . for ggh: 
# 



(A) 



no . 



27934000 e+02 
1 . 16637000 e-05 
1 . 17200000 e-01 
9 . 11876000 e + 01 
4 . 20000000 e + 00 
1.73300000e+02 
1 . 27500000 e+00 
MASS 

125.000000d0 
DISTRIB 



1 = light quarks at NLO , 
Standard Model inputs 
alpha_em - (-1) (MZ) SM MSbar 
G_Fermi 

SM MSbar 



2 = SM EW factor 



alpha_s (MZ) 
m_Z(pole) 
m_b (m_b ) 

# m_t(pole) 

# m_c (m_c ) 



# Higgs mass 



21 

22 
3 

31 
32 
4 

Block 
1 
2 

Block 
1 

Block 
1 
2 
3 
4 

Block 
1 
2 
3 

Block 
1 
2 
3 



30. dO 
100. dO 




distribution 



sigma_total , 1 
dsigma/dy, 3 
for pt and 
> ptmin , 2 
< ptmax 
in GeV 
in GeV 

rapidity - cut : = no , 1 = Abs [y] 
2 = Abs [y] > ymin , 3 = ymin < 
minimal rapidity ymin 
maximal rapidity ymax 
= rapidity , 1 = pseudorapidity 



pt -cut : 

minimal 
maximal 



(values 
= no , 1 = pt 
3 = ptmin < pt 
pt-value ptmin 
pt-value ptmax 



dsigma/dpt , 
d~2 sigma/ dy/dpt 
: 22 and 32) 
pt < ptmax , 



< ymax . 
Abs [y] 



ymax 



. 5d0 

1.5d0 # 
# 
SCALES 
1 . dO # 
1 . dO # 
REN0RMB0T # 
# 
PDFSPEC 

MSTW20081o68cl . LHgrid 
MSTW2008nlo68cl . LHgrid 
MSTW2008nnlo_asmzrange . LHgrid 



r enormal izat i on scale muR/mh 
factorization scale muF/mh 
Renormalization of the bottom 
m_b used for bottom Yukawa: 



= OS, 1 = MSbar(mb), 2 = MSbar(muR) 



name 
name 
name 



of 
of 
of 



pdf 
pdf 
pdf 



(lo) 

(nlo) 

(nnlo) 





VEGAS 
10000 
5 

10 

FACTORS 

0. dO 

1 . dO 
1 . dO 



# set number 



number of points 
number of iterations 
print : = no output 



1 = prettyprint , 10 = table 



factor for yukawa - coupl ings 

t 

b 
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Block SUSHI specifies the crucial input for SusHi, namely the model, the kind of Higgs 
boson to be considered (scalar or pseudo-scalajfl) , the type of collider, the center-of-mass 
energy, the perturbative order for gluon fusion and bottom-quark annihilation, and to which 
extent electro- weak corrections to gluon fusion should be taken into account. 

Block SMINPUTS contains the relevant SM input. We use the electro-magnetic coupling 
qeM) Fermi's constant Gp, and the Z-boson mass m z (entries 1,2,4) to calculate the W 
mass mw and the weak mixing angle sin 9w- The input value for a s (mz) given in entry 3 is 
used for renormalization-group (RG) running and RG transformations. We allow this value 
to be different from the one required by the PDFs which are specified further below. The 
latter is taken from LHAPDF and enters the calculation as the coupling parameter of the 
perturbative expansion of the cross section, for example Eq. (fT9"|) and Eqs. (|2~T|) . The charm- 
and bottom-quark masses (entries 8,5) are to be given in the MS scheme as m c (m c ) and 
trib{mb), while the top-quark mass (entry 6) is required in the on-shell scheme. In the SM, 
the Higgs mass is a free parameter and has to be provided in Block MASS, entry 1. 

Block DISTRIB controls cuts or distributions with respect to the transverse momentum 
Pt, the (pseudo-) rapidity y (tj), if desired. Note that differential cross sections (entry 1 
G {1,2, 3}) can only be obtained for gluon fusion; in this case, entries 22 and/or 32 specify 
the value of px and/or y (rj). Cuts are possible both for gluon fusion and bottom-quark 
annihilation; they are applied by setting entry 1 to ("total cross section"), and specifying 
entries 2,21,22 and/or 3,31,32. Note that entry 4 changes between rapidity y and pseudo- 
rapidity T). 

Block SCALES defines the renormalization and factorization scales relative to the Higgs 
mass. In accordance with Section 12.31 Block RENORMBOT offers different options for the 
renormalization of the bottom Yukawa coupling; three options are currently implemented: 
m\ e {mP s ,m^ s (m fe ),mf s (/i R )}. Block PDFSPEC contains the PDF sets in the notation 
of LHAPDF, consisting of the name of the PDF grid file, and the set number. Block VEGAS 
specifies integration parameters; note that distributions or cuts require higher numerical 
precision than the total cross section in order to reach comparable accuracy in the final 
result. Finally, Block FACTORS allows for additional factors in the Yukawa couplings of the 
fermions. We add that also charm-quark contributions can be taken into account by setting 
the corresponding factor to 1. Then the c-quark contributions at (N)LO are added using the 
on-shell value m^ s calculated from m c (m c ) as done for the on-shell bottom-quark mass. In 
case of the MSSM, for which a detailed prescription follows, the charm-quark contributions 
can be added as well. 



6 Apart from the scalar Higgs boson of the actual SM, SusHi also provides results for a pseudo-scalar 
Higgs-like particle whose coupling to fermions is obtained from the corresponding MSSM couplings by setting 
tan/3 = 1. 
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6.4-2. Minimal Supersymmetric Standard Model 

In case of the MSSM, the input file contains a number of additional Blocks. We show 
them here, together with the relevant modifications of the SM version: 



Block 
1 
2 

[ 

5 

Block 
3 

Block 

3 
11 
12 
23 
26 
43 
46 
49 
Block 

1 

2 

3 

4 

Block 
1 
2 
3 

Block 
1 
2 
3 

Block 
1 
2 
3 
4 
5 



SUSHI 
1 
1 
] 



2 

MINPAR 
5. dO 
EXTPAR 
800 . dO 
2006 . 66d0 
2006 . 66d0 
200 . dO 
300 . dO 
1000 . dO 
1000 . dO 
1000 . dO 
FEYNHIGGS 
0. dO 
200 . dO 
2006 . 66d0 
1000 . dO 
REN0RMB0T 

1 
1 

REN0RMSB0T 

2 





FACTORS 

0. dO # fac 

1 . dO # t 
1 . dO # b 
l.dO # st 
l.dO # sb 



# model : = SM , 1 = MSSM 

# = light Higgs (h) , 1 = pseudoscalar (A) , 2 = heavy Higgs (H) 

# order ggh : -1 = off, = LO , 1 = NLO , 2 = NNLO , 3 = " NNLO stop for h 

# SUSY breaking input parameters 

# tanb 



M_3 
A_t 
A_b 

mu in GeV 

M_A0 

M_Q3 

M_TR 

M_BR 

FeynHiggs specific input 

M_l 

M_2 

A in GeV (except for A_t , A_b) 
M_SUSY in GeV (except for M_Q3 , M_TR , 
Renormal izat i on of the bottom sector 
m_b used for bottom Yukawa: = OS , 
t an ( bet a ) - r es . of Y_b : = no. 



M_BR) 



= MSbar(m_b), 2 = MSbar(muR) 
1 = naive, 2 = full (for OS only) 



Delta_b : Take Delta_b from FeynHiggs: = no , 
# Renormal izat i on of the sbottom sector 



1 = 



# m_b : = OS . 

# A_b : = OS , 

# theta_b : = OS . 



1 = DRbar , 
1 = DRbar , 
1 = DRbar 



2 = dep 
2 = dep 



recommended : 
recommended : 
recommended : 



yes 

2 





tor for yukawa - couplings : c 



Entry 2 of Block SUSHI now distinguishes between the three MSSM Higgs bosons. En- 
try 5 allows to add approximated NNLO stop contributions for the light Higgs h. Block 
MINPAR, entry 3 defines the value of tan/3. Block EXTPAR fixes the parameters of the third 
family of quarks and squarks in the MSSM. If the Block FEYNHIGGS is present, SusHi has 
to be linked to FeynHiggs (see Section 16.3 j) which will then be used to calculate the Higgs 
masses from the parameters of that Block. In addition to the SM version, Block REN0RMB0T 
provides various ways of resumming t^-enhanced effects for the on-shell bottom Yukawa cou- 
pling, see Section [2731 For a running coupling (entry 1 e {1, 2}), resummation of those effects 
is always performed as shown in Eq. (I16fj) . Block REN0RMSB0T provides the choice between 
the various options of Tab.[TJ The alternative to Block FEYNHIGGS is the specification of 
the Higgs masses and the Higgs mixing angle a by hand; for example: 



Block ALPHA 

-2 . 58961078E-01 
Block MASS 

25 125 . 216431 E+00 

26 303 . 288802 E+00 
36 130 . 000000 E+00 



# mixing in Higgs sector 

# Higgs mass h 

# Higgs mass H 

# Pseudoscalar Higgs mass A 
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In this case, it is the user's responsibility to assure consistency of the Higgs mass and 



the other parameters. However, this option allows one to use H3 m |77f . for example, in order 
to take into account three-loop effects to the SUSY Higgs mass 76|~78j]. 

Three example input files can be found in the subfolder example, namely a SM input file 
and two MSSM input files, the latter two for the usage with and without FeynHiggs. 

6.4.3. Output file 

SusHi outputs the results of the calculation as well as some key parameters derived from 
the input in the same format as the input file. A typical example is shown here: 



Block SUSHIggh # Bon appetit 

1 1 . 65830627 E + 01 # ggh XS in pb 

Block SUSHIbbh # Bon appetit 

1 3.83211336E-01 # bbh XS in pb 
Block XSGGH # ggh MSSM-Cross sec. in pb 

2 1 . 21882475 E + 01 # NLO 
[ individual channels] 

Block XSGGHEFF # ggh MSSM -Cross sec. 

1 1.44925736E+01 # gghONLO SM 

2 1 . 82888376 E + 01 # gghONNLO SM 

3 5 . 78917317 E-02 # electroweak factor 
Block XSBBH # bbh MSSM -Cross sec. in pb 

1 5.30129833E-01 # L0 

2 4.79035691E-01 # NLO 

3 3.83211336E-01 # NNLO 

Block HGGSUSY # couplings of light Higgs h to 3. generation 
[ ] 

Block MASSOUT 

5 4.20000000E+00 # m_b ( m_b ) , MSbar 

25 1 . 25216431 E + 02 # MSSM -Mh in GeV 

[ SM masses / sbott om / st op masses] 

Block ALPHA # Effective Higgs mixing parameter 

-2 . 58961078E-01 # alpha 
Block STOPMIX # stop mixing matrix 
1 1 7.07918788E-01 # V_ll 

1 2 -7 . 06293841E-01 # V_12 

2 1 7 . 06293841 E-01 # V_21 

2 2 7.07918788E-01 # V_22 
[ ] 

Block AD 

3 3 2 . 00666000 E+03 # used Ab in GeV - def . accord, to scheme 
Block AU 

3 3 2 . 00666000 E + 03 # used At in GeV 



The main result for the gluon fusion cross section, containing all corrections specified by 
the user in the input file, is given as entry 1 in Block SUSHIggh, the one for bottom-quark 
annihilation in Block SUSHIbbh. Individual contributions to the cross sections are listed in 
Block XSGGH and Block XSBBH; their meaning should be obvious from the comment in the 
output file. Note that the results denoted "L0" etc. mean that the LO partonic cross section 
is convolved with the PDF set given in entry 1 of Block PDFSPEC in the input file. 

For gluon fusion, the Block XSGGHEFF contains the NNLO top-(s)quark results as ob- 
tained by ggh@nnlo, and the electro-weak correction factor as determined in Section l3~il 

In addition, Block HGGSUSY lists the non-resummed MSSM couplings of the quarks and 
squarks to the Higgs boson under consideration. Block MASSOUT gives the relevant SM and 
SUSY masses as well as the Higgs mass. Not shown above is the Block INTERNALMASSES, 
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which provides the different bottom masses entering the calculation of gluon-fusion cross 
sections, and SCALESOUT showing the renormalization and factorization scale as well as 
the value of q; s (/xr) taken from the PDF set at (N)LO. Finally all output files have the 
corresponding input file attached at the end. 



7. Conclusion 

In this article we described the Fortran code SusHi for the calculation of the cross section 
for Higgs production in gluon fusion and bottom-quark annihilation at hadron colliders. It 
works both in the SM and the MSSM, evaluates inclusive cross sections, distributions, and 
allows for kinematical cuts on the Higgs 4-momentum. It includes higher order QCD and 
electro-weak corrections and takes into account the effect from squarks and gluinos. 

SusHi allows one to choose among various renormalization schemes for the sbottom sector 
and the bottom Yukawa coupling, and includes the resummation of tan /3-enhanced effects. 
For the calculation of the Higgs mass in the MSSM SusHi can be linked to FeynHiggs. SusHi 
can be downloaded from Ref. 111. 
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Appendix A. Formulas: Higgs-squark couplings 

In this section we present the couplings of the three neutral Higgs bosons <fi of the MSSM 
to the quarks and squarks being implemented in SusHi. The relevant Feynman rules can be 
written in the form 




where v = 2mw/g — 1/v \f2Gp = + v^. The couplings of the Higgs boson <\> to 
the quarks q with respect to the SM Higgs boson coupling were already presented in Eq. ([I]) . 
The couplings g-^ of the squarks to the light and heavy Higgs can be split in the form 

9q,ij ~ 9q,ij + 9q,ij + 9q,ij ■ l A ^J 

In case of the light Higgs h we obtain for the couplings: 

h,EW _ EW 2 , EWJ h,EW _ EW 2 , JWJ /a q\ 

h,EW _ EW 2 , EW 2 ^,EW _ EW 2 , EW 2 /a 4 \ 



h.EW h,EW 1 / EW EW\ ft,EW h,EW L ( EW EW\ (\ r\ 

^,12 = %,21 =2^ ~ C 1 ) S2 ^ ^,12 =^,21 =2^ 2 ~° l ^ (A ' 5) 

%,n = -%, 22 = 9 %u = ~g~ K22 = --— *»s (A-6) 

%,12 = ^,21 = ^ C ^ V2 = V = ^ ^ (A ' 7) 

rffi = - (2 + 1 rfS - -- fa + ^T"^ 1 (A.8) 

y *> n S)3 \ 2m 2 2e V C/ ? I 2m 2 2e y v ; 

9 ^ ~ s p [ 2 2m\ S2e ^,22 " C/J ^ 2m 2 s 2, s J (A.9) 

%,12 " %,2! " - 2m? S2 % C2 ^ t - %,12 " %,21 " 2m 2 (A.10) 

Therein we made use of the abbreviations s x = sin x and c x = cos x and defined: 

^ = ~% (l - -n(a + /?) ofT - ^ (l - ^L) M« + ft (A-H) 

c ^ = -^f sin(a + ^ c ^ = if sin(a + 0) (A - 12) 

2.3 



The couplings to the heavy Higgs H are easy to obtain by the replacement a — > a — | in 
the previous formulas. In case of the CP-odd Higgs A the couplings are given by: 

^ii = ^ 22 = <ii = < 22 = ( A -!3) 

^-^^T S ^ + ^( 1 + ?) (A ' 14) 
< 12 = -< 21 = *^ ! ^r^^ + ^ C 1 + « (A. 15) 

We add that m& is partially interpreted as the bottom mass in the sbottom sector, namely 
where it is meant to be part of the Higgs-sbottom coupling. 
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